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' Abstract 

CN ■ 

The Kowalevski gyrostat in two constant fields is known as the unique example of 
an integrable Hamiltonian system with three degrees of freedom not reducible to a fam- 
ily of systems in fewer dimensions and still having the clear mechanical interpretation. 
The practical explicit integration of this system can hardly be obtained by the existing 
techniques. Then the challenging problem becomes to fulfil the qualitative investigation 
based on the study of the Liouville foliation of the phase space. As the first approach to 
topological analysis of this system we find the stratified critical set of the momentum map; 
this set consists of the trajectories with number of frequencies less than three. We obtain 
the equations of the bifurcation diagram in three-dimensional space. These equations have 
the form convenient for the classification of the bifurcation sets induced on 5-dimensional 
■ iso-energetic levels. 
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[Z] ; 1 Introduction 

m ' 

The famous integrable case of S. Kowalevski of the motion of a heavy rigid body about a fixed 
point pQ has received several generalizations. Some of them suppose restrictions to submanifolds 
in the phase space (partial cases), others are far from mechanics involving potential functions 
on the configuration space SO(3) with singularities. The most essential generalization having 
the clear mechanical sense was found by A.G. Reyman and M.A. Semenov-Tian-Shansky in the 
work 2j. The authors introduce the dynamical system on the dual space of the Lie algebra 
e(p, q) of the Lie group defined as the semi-direct product of SO(p) and q copies of K p . Such 
systems are known as the Euler equations on Lie (co)algebras [3J. The case p — 3, q = 2 
corresponds to the Euler-Poisson equations of the motion of a gyrostat in two constant fields. 

For a rigid body without gyrostatic momentum, the model of two constant fields was intro- 
duced by O.I. Bogoyavlensky [3J. The physical object can be either a heavy electrically charged 
rigid body rotating in gravitational and constant electric fields, or a heavy magnet rotating in 
gravitational and constant magnetic fields. The corresponding equations are Hamiltonian on 
the orbit of coadjoint action on e(3, 2)* of the Lie group defined as the semi-direct product 
SO(3) x (M 3 (g) R 3 ). The typical orbit is diffeomorphic to TSO{3) = R 3 x SO(3). Therefore, 
the gyrostat in two constant fields is the Hamiltonian system with three degrees of freedom. 
Bogoyavlensky 3] suggested the conditions of the Kowalevski type and found the analogue 
of the Kowalevski integral K for the top in two constant fields. H. Yehia [1] generalized this 
integral for the Kowalevski gyrostat in two constant fields. Almost simultaneously with Yehia, 
I.V. Komarov [5] and L.N. Gavrilov [5] proved the Liouville integrability of the Kowalevski gyro- 
stat in the gravity field. But for two constant fields the Kowalevski gyrostat was not considered 
integrable due to the fact that the existence of the second field destroys the axial symmetry 
of the potential and, consequently, the corresponding cyclic integral. Finally, Reyman and 
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Semenov-Tian-Shansky 2 found the Lax representation with a spectral parameter for the fam- 
ily of Euler equations on e(p, q)* . For e(3, 2) this representation immediately gave rise to the 
new integral for the Kowalevski gyrostat in two constant fields. For the classical Kowalevski 
top this integral turns into the square of the cyclic integral. 

The Kowalevski gyrostat in two constant fields does not have any explicit symmetry groups 
and, therefore, is not reducible, in a standard way, to a family of systems with two degrees of 
freedom. Phase topology of such systems has not been studied yet. The theory of n-dimensional 
integrable systems started in [7J is not illustrated by an application to any real irreducible 
physical or geometrical problem with n > 2. 

In the paper [8] , the authors give a detailed exposition of the results of [2] as well as a study 
of the algebraic geometry of the Lax pair for the generalized Kowalevski system. They announce 
the possibility of its integration by the finite-band techniques and fulfil such integration for the 
classical top. For two constant fields the integration of the Kowalevski top is not given up-to- 
date. The problem of the Kowalevski gyrostat motion in two constant fields is not studied at all. 
The technical difficulties here are extremely high. It is not likely that, in the general regular case, 
the analytical solutions can be obtained having the form useful for the qualitative topological 
analysis or the computer simulation. However, there is a good experience of studying the critical 
subsystems, i.e., the systems with two degrees of freedom induced on 4-dimensional invariant 
submanifolds of the phase space. For the Kowalevski top in two constant fields we have now the 
complete description of all singularities of the momentum map [9] , [10] , [IT] , [12] , [13] , [14] and 
the classification of the bifurcation diagrams for the restriction of this map to 5-dimensional 
iso-energetic surfaces [15], [16], [17] . This result is a necessary and highly complicated part 
of the study of Liouville foliation of the integrable system and shows the actual need in the 
generalization of the Liouville invariants theory [15] for the dimensions greater than two. 

The present paper contains similar results for the Kowalevski gyrostat in two constant 
fields. The 6-dimensional phase space is stratified by the rank of the momentum map. We find 
the equations of invariant submanifolds on which the induced systems are Hamiltonian with less 
than three degrees of freedom (critical manifolds of rank 0,1, or 2). We straightforwardly prove 
that the image of these critical manifolds (the bifurcation diagram) lies in the discriminant set of 
the algebraic curve of the Lax representation given in [8]. Moreover, the spectral parameter on 
the Lax curve is explicitly expressed in terms of the constant s of the additional partial integral 
arising on the critical submanifolds. It then follows that the equations of the surfaces containing 
the bifurcation diagram are written in the parametric form such that the parameters are the 
energy constant h and the constant s of the partial integral. Fixing the value of h we come to 
explicit equations of the bifurcation diagrams induced on iso-energetic levels. The problem of 
classification of these diagrams seems quite complicated due to the existence of several physical 
parameters. Nevertheless, it is certainly solvable with the help of contemporary computer 
programs of analytical calculations. 

First we show that the number of physical parameters for the gyrostat in two constant fields 
can be reduced by a simple procedure, which may be called the orthogonalization of the fields. 
More precisely, for the problems of gyrostat motion there exists a group of diffcomorphisms of 
the phase spaces (mentioned above orbits of the coadjoint action) that is an equivalence group 
for the corresponding dynamical systems. It appears that each equivalence class contains a 
problem with an orthonormal pair of radius vectors of the centers of forces application and with 
an orthogonal pair of the intensity vectors. Such force field is characterized by only one essential 
parameter — the ratio of the modules of the intensity vectors. For a dynamically symmetric 
gyrostat having the centers of forces application in the equatorial plane, the orthogonalization 
procedure along with the appropriate choice of the measure units leave, in addition to the 
forces ratio, only two physical parameters of the body itself, namely, the ratio of the equatorial 
and axial inertia moments and the non-zero axial component of the gyrostatic momentum. In 
the generalized Kowalevski case the first of them equals 2. Thus, the whole problem has, in 
fact, two essential parameters. In particular, each of the critical four-dimensional submanifolds 
found below provides a two-parametric family of completely integrable Hamiltonian systems 
with two degrees of freedom. 
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2 Gyrostat equations and parametrical reduction 



Consider a rigid body B rotating around a fixed point O. Choose a trihedral at O moving along 
with the body and refer to it all vector and tensor objects. Denote by eie2e3 the canonical unit 
basis in R 3 ; then the moving trihedral itself is represented as Oeie2e3. Let to be the vector 
of the angular velocity of B. Suppose that B is bearing an axially symmetric rigid rotor B' 
rotating freely around its symmetry axis fixed in B. Such system of two bodies is the simplest 
model of a gyrostat. The notion of a gyrostat was introduced by N.E. Zhukovsky [19] for a body 
having cavities totally filled with homogeneous fluid. Both models have the common feature 
usually taken as the definition of a gyrostat: the total angular momentum of such system is 
M = 111) + A, where the inertia tensor I and the vector A (called the gyrostatic momentum) are 
constant with respect to the moving trihedral. Using the term "gyrostat" we always suppose 
A ^ 0. In the case A = we use the terms "rigid body" or "top" instead. The top is usually 
supposed to have a dynamical symmetry axis. 

Let denote the moment of external forces with respect to O (the rotating moment). 
Constant field is a force field inducing the rotating moment of the form r x a with constant 
vector r and with a corresponding to some physical vector fixed in inertial space; r points from 
O to the center of application of the field, ol is the field intensity. 

For two constant fields the rotating moment is Mi? = ri x a + r2 x (3 with ri , T2 constant 
in the body and a, (3 corresponding to the vectors fixed in inertial space. Obviously, Mp can 
be represented as the moment of one constant field if either ri x r*2 = or a x f3 = 0. Suppose 
that 



Two constant fields satisfying (|2.ip are said to be independent. 

The equations defining the respective evolution of M,a,/3 in two constant fields are 



These equations are Euler equations in the space R 9 (M, a, f3) considered as the dual space 
to the semi-direct sum so(3) + (R 3 <£> R 3 ). The Lie-Poisson bracket applied to the coordinate 
functions yields 



Such bracket is non-degenerate on each orbit of the coadjoint action. The orbits are defined by 
the geometric integrals (common level of the Casimir functions) 



If en > 0, c 2 2 > 0, c\ 2 < C11C22, then the orbit in R 9 is diffeomorphic to R 3 x 50(3), and the 
induced Hamiltonian system has three degrees of freedom (see [3J, [H] for the details). From 
physical point of view the constants cn,C22,ci2 characterize the force fields intensities. Along 
with the coordinates of rx,T2 in the moving frame, we have 9 parameters of the interaction of 
the body with the external forces. We now show how to reduce this number. 
Introduce some notation. 

Let L(n, k) be the space of n x fc-matrices. Put L(k) = L(k, k). 

Identify R 6 = R 3 x R 3 with L(3, 2) by the isomorphism j that joins two columns 



rj x r 2 / 0, ax/3/0. 



(2.1) 




(2.2) 



{Mi, Mj} = EijkMk, {Mi, aj} = eijkak, {Mi, ft} = e ijk (3k 
{a t , a 3 ) = 0, {a. t , = 0, {ft, ft} = 0. 



(2.3) 



a • ct = en, /3-/3 = c 2 2, a-/3 — ci 2 . 



A = j'(ai,a 2 ) = ||ai a 2 || G L(3, 2), 



ai,a 2 G R 3 . 



For the inverse map, we write 



j- 1 (A) = (c 1 (A),c 2 (A))gR 3 xR 3 



A G L(3,2). 
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If A, B £ L(3, 2), a £ R 3 , by definition, put 

2 

Ax B =f^ Ci{A) x Cl (B) £ M 3 ; 

axi = j(ax ci(A),a x c 2 (A)) e £(3,2). 

Lemma 1. let A e 50(3), D e GL(2, R), a e R 3 , A, B e 1(3, 2). Then 

A(A xB) = (AA) x (AS); {AD- 1 ) x [BD T ) =AxB; 
A(a x A) = (Aa) x (AA); a x (AD) = (a x A)D. 

The proof is by direct calculation. 

In notation (|2.4|) we write Eqs. (|2.2p in the form 



(2.4) 



1^ = {lu + \) x u + A x U, j = -uxK (2.5) 

Cit tit 

Here A = j( r ij r 2) is a constant matrix, £/ = j(ct,/3). The phase space of (|2.5p is {(<*>, t/)} = 
R 3 x L(3,2). 

In fact, [/ in (12. 5p is restricted by the geometric integrals; i.e., for some constant symmetric 
matrix C £ L(2) 

U T U = C. (2.6) 

Let O be the set defined by Eq. (|2.6[) in Z/(3, 2). In order to emphasize the O-dependence, we 
write O = 0(C). 

Let <P = (I, A, A, C) denote the complete set of constant parameters of the problem. Denote 
by Xtp the vector field on R 3 x 0(C) induced by (12. 5| . Given the set ^3, the problem of motion 
of the gyrostat in two constant fields described by the dynamical system ATtp will be called, for 
short, the DG-problem. 

Associate to A g 50(3), D £ GL(2,R) the linear automorphisms W(A,D) and ip(A,D) of 
R 3 x L(3, 2) and L(3) x M 3 x L(3, 2) x L{2) 



¥(A, D)(u>, £/) = (Aw, AC/Z? T ), 

^(A, £>)(!, A, A, C) = (AIA T , AA, AAD -1 , DCD T ). 



(2.7) 



Eqs. (EU) and (g7F|) imply *(A,L>)(JR 3 x 0(C)) = R 3 x 0(DCD T ). Using Lemma Q] we 
obtain the following statement. 

Lemma 2. For eac/i (A, 13) e 50(3) x GL(2,R), we have 

9(A,D).(Xy{v))=X^ D)m (9(A,D){v)), v £ R 3 x 0(C). 

Thus, any two DG-problems determined by the sets of parameters and ip(A, D)^) are 
completely equivalent. 

Let us call a DG-problem canonical if the centers of application of forces lie on the first 
two axes of the moving trihedral at unit distance from O and the intensities of the forces are 
orthogonal to each other. 

Theorem 1. For each DG-problem with independent forces there exists an equivalent canonical 
problem. Moreover, in both equivalent problems the centers of application of forces belong to the 
same plane in the body containing the fixed point. 

Proof. Let the DG-problem with the set of parameters = (I, A, A, C) satisfy ()2.1|) . This 
means that the symmetric matrices A* = (A T A)- 1 and C are positively definite. According to 
the well-known fact from linear algebra, A* and C can be reduced, respectively, to the identity 
matrix and to a diagonal matrix via the same conjugation operator 

DA*D T = E, DCD T = diag{a 2 , 6 2 }, D £ GL(2, R), a, b £ R+. 
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Then c\(AD~ 1 ) and c 2 (AZ?~ 1 ) form an orthonormal pair in R 3 . There exists A E SO{2>) such 
that Aci(^4D -1 ) = e, (i = 1, 2). The first statement is obtained by applying Lemma [5] with the 
previously chosen A, D to the initial vector field X<p. 

To finish the proof, notice that the transformation A i— > AD^ 1 preserves the span of C\(A), 
c 2 (A). The matrix A in (|2.7[) stands for the change of the moving trihedral. Therefore, if a € R 3 
represents some physical vector in the initial problem, then Aa is the same vector with respect 
to the body in the equivalent problem. □ 

Remark 1. The fact that any DG-problem can be reduced to the problem with one of the 
pairs i"i,r2 or c*,/3 orthonormal is obvious. Simultaneous orthogonalization of both pairs was 
first established in for a rigid body and crucially simplifies all calculations. 

It follows from Theorem [1] that, without loss of generality, for independent forces we may 
suppose 

ri = ei, r 2 = e 2 , (2.8) 
a-a = a 2 , f3-f3 = b 2 , a ■ [3 = 0. (2.9) 

Change, if necessary, the order of ei,e2 (with simultaneous change of the direction of e^) 
to obtain a ^ b > 0. 

Consider a dynamically symmetric top in two constant fields with the centers of application 
of forces in the equatorial plane of its inertia ellipsoid. Choose a moving trihedral such that 
Oe^ is the symmetry axis. Then the inertia tensor I becomes diagonal. Let a = b. For any 

S SO(2) denote by 6 € SO(3) the corresponding rotation of K 3 about Oe 3 . Take in (j2~7f 
A = 0, D = 0. Under the conditions (|2.8|) . (|2.9|) . ip — Id and iff becomes the symmetry group. 
The system (|2.5j) has the cyclic integral Iu> ■ (a 2 e3 — a. x /3). Therefore it is possible to reduce 
such a DG-problem to a family of systems with two degrees of freedom. For the analogue of 
the Kowalevski case this system becomes integrable [I]. 

Let us call a DG-problem irreducible if, in its canonical representation, 

a>6>0. (2.10) 

The following statements are needed in the future; they also reveal some features of a wide 
class of DG-problems. 

Lemma 3. In an irreducible DG-problem, the body has exactly four equilibria. 

Proof. The set of singular points of (|2.5|) is defined by lj — 0, A x U — 0. For the equivalent 
canonical problem with (|2.8p we have 

ei x a + e 2 x [3 = 0. 

Then the four vectors ei, a, e 2 , (3 are parallel to the same plane and |ei x a\ = |e 2 x (3\. Given 
()2.10|) . this equality yields 

a = ±aei, (3 = ±be 2 . (2.11) 

Thus, in the canonical irreducible system, an equilibrium takes place only if the radius vectors 
of the centers of application are parallel to the corresponding fields intensities. □ 

Note that the existence of the gyrostatic momentum does not change the equilibria. There- 
fore, the result here is the same as in the case of a rigid body in two constant fields [15] . 

Lemma 4. Let an irreducible DG-problem in its canonical form have the diagonal inertia tensor 

1 = diag{7i, / 2 , J3} and X = 0. Then the body has the following families of periodic motions of 
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pendulum type 



p f w = ip'ei, a = ±aei, [3 = b(e 2 cos ip - e 3 sin^), (2 12) 

1 ' 1 Ii<f" = —bsinip; ^ ' ' 

p f w = <^'e 2 , /3 = ±6e 2 , a = a(ex cos tp + e 3 sin ip), (2 13) 

2 ' [ I 2 tp" = -asin<^; 

{u> = </?'e3, a x j3 = ±a6e3, 
a = a(ei cosy; — e 2 siny), (3 = ±6(ei sin </? + e 2 cos </?) , (2-14) 
73<£>" = — (a ± 6) sin 

If X but X = Xei for some i — 1,2,3, then the only family remained is Pi with the 
corresponding index. 

The proof is obvious. The families (|2.12[) - (12.14[) were first found in [TT] (the case A = 0). 
Note that for two constant fields these families are the only motions with a fixed direction of 
the angular velocity. In particular, the body in two independent constant fields does not have 
any uniform rotations. 



3 Critical set of the Kowalevski gyrostat 

Suppose that the irreducible DG-problem has the diagonal inertia tensor with the principal 
moments of inertia satisfying the ratio 2:2:1, the gyrostatic momentum is directed along the 
dynamical symmetry axis A = Ae 3 and the centers of the fields application lie in the equatorial 
plane ri_Le3, r 2 _Le3. These are the conditions of the integrable case [2] of the Kowalevski 
gyrostat in two constant fields. The orthogonalization procedure in this case does not change 
the e3-axis and we obtain (|2.8[) . (|2.9p . Choosing the appropriate units of measurement, represent 
Eqs. (|2.5|) in the form 



2wi = cj 2 (w 3 - A) + /3 3 , 2cj 2 = -wi(w3 - A) - a 3 , lj 3 = a 2 - (3i 
di = a 2 uj 3 - a 3 u> 2 , $i = - 
a 2 = a^uji - aiw 3 , $2 = (hui - 0iW3, 
a 3 = a X L)2 - ce 2 u>i, $3 = (3\U} 2 - @2U>i. 



(3.15) 



The phase space is P 6 = R 3 x O, where O c I 3 x t 3 is defined by (|2.9p : O is diffeomorphic to 
50(3). 

The complete set of the first integrals in involution on P e includes the energy integral H , 
generalized Kowalevski integral K [3] , [4] , and the integral G found in 2J . After the parametrical 
reduction, these integrals are 

H = uj 2 + u)\ + -W3 - «i - [32, 

K = (wf - uj\ + ai - (3 2 ) 2 + (2wiw 2 + a 2 + (3 1 ) 2 + 

+ 2X[(uj3 - X)(oj 2 +uj 2 2 ) + 2uj ia3 + Zwalh], 
G = + M 2 ) + i(«3 - A)M 7 - b 2 ai - a 2 l3 2 . 

Here M a = (Iw + X)-a, Mp = (Iu> + X)-f3, Af 7 = (Iw + A) -(a x 0). 
Introduce the momentum map 

J = G x K x H : P 6 ^R 3 (3.16) 

and denote by € C P 6 the set of critical points of J. By definition, the bifurcation diagram of 
J is the set Eel 3 over which J fails to be locally trivial; S defines the cases when the integral 
manifolds 

J c = J- 1 {c) 1 c=(g,k,h)eR 3 
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change its topological (and smooth) type. To find £ and E is the necessary part of the global 
topological analysis of the problem. 

It follows from Liouville-Arnold theorem that for c (ji E the manifold J c , if not empty is the 
union of three-dimensional tori. The considered Hamiltonian system on P 6 is non-degenerate at 
least for small enough values of b. Therefore the trajectories on such tori are almost everywhere 
quasi-periodic with three independent frequencies. The critical set £ is preserved by the phase 
flow and consists of the trajectories having less than three frequencies. We call these trajectories 
the critical motions. The set £ is stratified by the rank of J. Let £.,■ = € £ : rank J(C) = j} 
(j = 0, 1, 2). It is natural to expect that £j consists of the Liouville tori of dimension j and the 
image J(£j), as a subset of E, is a smooth surface Ej of dimension j. More precisely, for each 
j ^ 2 we have to take 

3-1 
i=0 

Then, as a whole, we may consider E as a two-dimensional cell complex, Ej as its j-skelcton. 
For j = 1, 2 we will have 9Ej C Ej_i. 

For c € E2 the set J c n £ consists of two-dimensional tori. Take the union of such tori 
over the values c from some open subset in E2. The dynamical system induced on this union 
will be Hamiltonian with two degrees of freedom. Vice versa, let M be a submanifold in P 6 , 
dim M — 4, and suppose that the induced system on M is Hamiltonian. Then obviously M C £. 
This speculation gives a useful tool to find out whether a common level of functions consists of 
critical points of J. 

Lemma 5. Consider a system of equations 

fi =0,...,/ 2fc = (3.17) 

on a domain W open in the phase space P 2n of the integrable Hamiltonian system X. Let 
M CW be the set defined by ]3.17[ . Suppose 

(i) fi, . . . , fik are smooth functions independent on M; 

(ii) Xf x = 0, . . . , Xf 2k = on M; 

(Hi) the matrix of the Poisson brackets is non- degenerate almost everywhere on 

M. 

Then M consists of critical points of the momentum map. 

Proof. Conditions (i), (ii) imply that M is a smooth (2n — 2fc)-dimensional manifold invariant 
under the restriction of the phase flow to the open set W. Condition (iii) means that the closed 
2-form induced on M by the symplectic structure on P 2n is almost everywhere non-degenerate. 
Thus the flow on M is almost everywhere Hamiltonian with n— k degrees of freedom. It inherits 
the property of complete integrability. Then almost all its integral manifolds consist of (n — k)- 
dimensional tori and therefore lie in the critical set of the momentum map. Since M is closed 
in W and the critical set is closed in P 2 ™ , we conclude that M totally consists of the critical 
points of the momentum map. □ 

Remark 2. In our case n = 3 and the above lemma is applied in the situations when k = 1 
or k — 2. The critical set and the bifurcation diagram of the map H3.16[) in the case A = are 
known. The critical set is described by one system of the type p,17[) with k = 2 and three 
systems of the type Q3.17P with k = 1. The complete presentation of these results and the list of 
publications are given in [TJ] , [T7] . Except for the partial integrable case of Bogoyavlensky [3] 
(case K = 0), all of the critical subsystems have been either explicitly integrated or reduced to 
separated systems of equations [T3], [H], [TO] . 

Introduce the change of variables [TU] based on the change given by S. Kowalevski and on 
the Lax representation [2] (i 2 = — 1) 

xi = (ai - f3 2 ) + i(a 2 + x 2 = (a\ - f3 2 ) - i{&2 + Pi), 

yi = (ai+f3 2 )+i(a 2 -f3i), y 2 = {a\ + f3 2 ) - i(a 2 - fa), (3 18) 

zi = «3 + ifa, Z2 — «3 - ifa, 
W\ = U\ + iuj 2 , w 2 = wi — iuj 2 , W3 = cj 3 . 
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Then Eqs. lETT5jl yield 



2w[ = -wi(w 3 - A) - zi, 2w' 2 = w 2 (w 3 - A) + z 2 , 2w 3 = y 2 -y 1: 
x\ = -X1W3 + Z1W1, x' 2 — x 2 w 3 — z 2 w 2 , 
y[ = ~yiw 3 + z 2 wi, y 2 = y 2 w 3 - ziw 2 , 
2z[ = x\u>2 — V2W\, 2z 2 = —x 2 w\ + y\w 2 . 



(3.19) 



Here prime stands for d/d(it). 

Consider (|3.18[) as the map M. 9 — ► C 9 and denote its image by V 9 . Eqs. (|2.9p of the phase 
space P 6 in V 9 take the form 

zf+x iy2 =r 2 , z^ + x 2yi =r 2 , (3.20) 
xix 2 + yiy 2 + 2z\z 2 = 2p 2 . (3-21) 

Here we introduce the positive constants 

p = \/ a 2 + b 2 , r = y/ a 2 — b 2 . 

Using (|3.20[) and p.2ip . express the first integrals in new coordinates, 

1 2 1. 

H = wiw 2 + -w 3 - -{yi + 2/2), 

K = (w 1 + xi)(w 2 + x 2 ) + 2X(wiw 2 w 3 + z 2 wi + zxw 2 ) — 2X 2 wiw 2 , 

G = j(p 2 - XiX2)wl + ^-{x 2 z 1 wi + x 1 z 2 w 2 )w 3 + 

1 1 (3-22) 

+ -(x 2 W! + yiw 2 )(y 2 W! + x x w 2 ) - -^P 2 (yi + 2/2) + 

+ jr 2 (xi + x 2 ) + ]^\{ziz 2 w 3 + y 2 z 2 wi + yiZiw 2 )+ 
+ \^ 2 {P 2 ~ 2/12/2)- 

Let / be an arbitrary function on V 9 . For brevity, the term "critical point of /" will always 
mean a critical point of the restriction of / to P 6 . Similarly, df means the restriction of the 
differential of / to the set of vectors tangent to P 6 . While calculating critical points of various 
functions, it is convenient to avoid introducing Lagrange's multipliers for the restrictions (|3 . 20[) 
and (|3~2T) . 

Lemma 6. Critical points of a function f on V 9 , in the above sense, are defined by the system 
of equations 

(3.23) 

where 

Xi 
Xi 

x 5 
x 6 

Indeed, six vector fields Xi are tangent to P 6 and linearly independent at any point of P 6 . 
The following two propositions define the strata £0 and £1 of the critical set. 

Proposition 1. The set £q consists exactly of the four equilibria existing in this problem. 



Xif 


= (i - 


= 1,...6), 
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Proof. The condition of zero rank of the momentum map at a point £ £ P e supposes, in 
particular, that dH — 0. Then £ is the point of equilibrium and it follows from Lemma [3] that 
C is one of the points (|2.1ip . Using the complex variables we have 

wi = w 2 = w 3 = 0, z\ = z 2 = 0, 
xi = x 2 = eia - e 2 b, yi = y 2 = E\a + e 2 b (ei = ±1, e 2 = ±1). 

Use Eqs. flXH]) with / = K and / = G to obtain that dK(() = and dG(() = 0. Therefore, 
rankJ(C)=0. □ 

Note that in classical problems of the rigid body dynamics with an axially symmetric force 
field, the rank of the momentum map is everywhere not less than 1 due to the regularity of the 
cyclic integral. In our case, all equilibria are non- degenerate (in the Morse sense) critical points 
of the Hamilton function (see [15)). Therefore, these points are critical for any first integral of 
the system. 

It is essential that in the sequel A ^ 0. 

Proposition 2. The set <L\ is completely defined by the condition 

reink{dK,dH} = 1 

and consists of the points of the following periodic trajectories: 

1) pendulum motions \2. 11$ ; 

2) motions defined by the equations 

wi = q(w)y/w, w 2 = -7^-, w 3 = -w, (3.24) 
q(w) a 

X X = — [r 2 X 2 a 2 -(X 2 +a)uq 2 (w)w], 
cni 

lr2\22 /\2, \ w i 

x 2 = — [ r A a — (A + a) u 



o~u 1 ' q 2 (w) ' 

a r 4 X 2 a. r 2 X 2 . 
yi = cr(l + — — ) H q (w)w, 



2/2 = 



A 2 M' ' u 

a r i X 2 a, r 2 X 2 w 



A 2 u 2 u <7 2 (i«)' 

r 2 Xa Jw X 2 + a _ 

Zi = -— H — q(w)Vw, 

u q(w) X 

r 2 Xa !— X 2 + a Jul 

z 2 = q(w)Jw H — . 

u X q(w) 

Here q(w) is the root of the equation q 4 — 2Q(w)q 2 + 1 = 0, where 

. an 3 + (A 2 + a)[X 2 w 2 + a 2 (2w - a)]u 2 + r 4 AV 4 

Q(w) = 2 2 2I 2 — : ; 3.26 

2r z X z a z (X z + a)uw 

a, u are constants satisfying the equation 

A 2 (A 2 + a) 2 u 5 + (A 2 + <j)[2p 2 X 4 - (A 2 + afa]cju A + . . 

+r 4 X 6 a 2 u 3 + 2r 4 X 4 a 4 (X 2 + a) 2 u 2 - r 8 X 8 a 6 = 0. [ - 6M > 

The evolution w(t) is defined by the equation 

(|) J = -^P + (»)P_H (3-28) 



where 



D ., 2 n 2 u±r 2 X 2 a[u 3 -(X 2 +a)a 2 u 2 +r 4 X 4 a 3 } , on , 

P ± (w)=w 2 + 2a 2 — w + ^ ' L. 3.29 

X z u [X z + a)X z u z 
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Proof. It follows from above that dH ^ at the points of (£1. Then to investigate the depen- 
dence of the functions K and H it is sufficient to introduce the function with one Lagrange's 
multiplier a. Write Eqs. P~2"3"j) with / = K - 2aH, 



(w\ + xi)w 2 + A[zi + wi(w 3 - A)] - awi = 0, , , 

(u>§ + x 2 )w\ + X[z 2 + w 2 (w3 - A)] - o~w 2 = 0, 

\w\w 2 — o~W3 = 0, (3.31) 

(w\ + x\)z 2 - \{w 2 xi + wiyi) + azi = 0, , , 

(w\ + x 2 )zi - \{w\x 2 + w 2 y 2 ) + az 2 = 0, 

x\w\ — x 2 w\ + a(y\ — y 2 ) = 0. (3.33) 

First consider the critical points of the function K. For this purpose we must put a = 0. 
Eq. (j33T|) gives Wi = w 2 = 0. Then Eqs. ([330")) imply z x = z 2 = 0. Eqs. (j33^) and ([3331) 
become identities. The same values satisfy Eqs. p.23[) if we take / = 4G + (xix 2 — y\y 2 )H. 
Therefore, dK = and 4dG + (xix 2 — y\y 2 )dH — 0. Since dH ^ 0, it means that rank J = 1. 
The initial variables on the corresponding trajectories are u)% = lo 2 = 0, 0:3 = $3 = 0. Substitute 
these values to Eqs. f|3. 15[) to obtain the solutions (|2. 14|) . 

Let a 7^ 0. The equilibria of the system are already excluded. Then it follows from (|3.31|) 
that wiw 2 7^ 0. Satisfying (|3.31j) . introduce new variables w,q as shown in (I3.24| . Four 
equations (|3.30|) . p.32p form the linear system in yi,yi, Zi, Z2, from which we obtain these 
variables as the functions of x\,x 2 ,w,q identically satisfying (|3.33|) . Denote 

u= (w-a) 2 {\ 2 +a)-ax lX2 . (3.34) 

Then Eqs. (|3.20p are easily solved for xi,x 2 as the functions of w,q,u. As a result we obtain 
the expressions (|3.25[) . Let 

<H<' 2 + ?»- 

Then the substitution of Xi, x 2 from (|3.25p back to (|3.34p gives (|3.26p . The last unused equation 
(|3.2ip provides the relation (|3.27p between u and the constants A, a. It shows that u defined 
as (I3.34j) appears to be a constant. 

Thus, all phase variables are expressed via one variable w, for which from (|3.19p we find 
the differential equation (|3.28|) . Note that due to (|3.29p the solutions are elliptic functions of 
time. 

To finish the proof, we need to show that at the points of the trajectories found we really 
have rank J = 1, i.e., the linear dependence of dK and dH implies the linear dependence of dG 
and dH. Indeed, Eqs. (PQgj) with 

f = 2G-{p 2 + ^±^u)H 
Y fl- 
are satisfied both by (j2~T4")l and by (|3~2"4"]l . PT23)) . Therefore, r&nk{dG, dH} = 1 and, conse- 
quently, rank{ dK, dG, dH} = 1. □ 

The following statement describes one of the critical subsystems in £2- 

Proposition 3. The system \3. 19\i has the four- dimensional invariant submanifold defined 
by the equations 

Ui =0, U 2 = 0, (3.35) 



where 

TT y2W\ + xiw 2 + zi(w 3 + A) x 2 w± + y±w 2 + z 2 (w 3 + A) 

Ui — , 

Wi w 2 

U 2 = w 1 w 2 U[. 

The Poisson bracket {Ui,U 2 } is non-zero almost everywhere on this submanifold. 



(3.36) 
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Proof. The derivative U 2 in virtue of (|3.19j) is proportional to U\, i.e., (I3.35|) implies U 2 = 0. 
Therefore the set p. 351) is invariant. 
Consider the function 



-4 



y 2 wi + xiw 2 + zi(w 3 + A) x 2 wi + y\w 2 + z 2 (w 3 + A) 



On D* we obtain 



W\ w 2 



, w 1 z 2 + w 2 z 1 +w 1 w 2 {wz~ X) TT _ n 
b — U\ — U. 

OWiW 2 



Therefore, S is a partial integral of the induced system. Eliminate 2/1,2/2 with the help of 
Eqs. (|3.35[) and present S in a more simple form 

s _ x 2 z 1 w 1 + xiz 2 w 2 + z 1 z 2 (w 3 + A) 
2wiw 2 (w 3 - A) 



Now the Poisson bracket of U\ and U 2 calculated under the rules defined by (|2.3p is expressed 
in terms of the energy constant h and the constant s of the integral (|3.37[) in the following way 



{Ui,U 2 } 



4 



•s 



3s 4 -2s 3 (fr- — ) + 



X\ p 4 -r 4 ^ 



2' 4 



Obviously, the right part of it is a ratio of polynomials not identically zero onD,. Therefore the 
set {Ui, U2} = has codimension 1 in D*. In particular this set is of zero measure in £)*. □ 

Remark 3. If A = 0, then the manifold D* turns into the phase space of the Hamiltonian 
system with two degrees of freedom studied in [14 . The geometrical characteristic of the 
motions in this system is the condition 

M-« M-/3 

— — = — — = const, 
M-ei M-e 2 

where M = la; is the angular momentum vector. The system (|3.35|) . (|3.36p is found from the 
same condition given that here M = Ilj + A. 

The following theorem completes the description of the critical set of the momentum map 
for the gyrostat. 

Theorem 2. The set of critical points of the momentum map 0(3.16} consists of the following 
subsets in P 6 : 

1 ) the set £ defined by the system 

wi =0, w 2 = 0, z 1 =0, z 2 = 0; (3.38) 

2) the set 9t defined by the system 

Ft =0, F 2 = 0, (3.39) 

where 

F\ = (wiw 2 + Xw 3 )(w 2 x 1 + Azi)Ayi — 

-w 2 (w\ + Xl)(x 2 ZlWl + XlZ 2 W 2 - XlX 2 W 3 + 2ziz 2 \)- 

-x 2 (wiw 3 + zi)(wizi - xiw 3 )X + (xiwl ~ 2ziw\w 3 - zf)z 2 X 2 , 
F2 = (wrw 2 + Xw 3 )(wix 2 + Xz 2 )Xy 2 - 

-Wi(»2 + X 2 ){x 2 ZiW\ + X\Z 2 W 2 — X\X 2 W 3 + 2zxz 2 X) — 

-xi(w 2 w 3 + z 2 )(w 2 z 2 - x 2 w 3 )X + (x 2 wl - 2z 2 w 2 w 3 - z 2 )ziX 2 ; 

3) the set O defined by the system 

Ri =0, R 2 = 0, (3.40) 
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where 

R\ = [j/iW 2 + X2W1 + z 2 (w 3 + \)]wi(w3 — A) + 

+ x 2 z\W\ + xiz 2 w 2 + ziz 2 {wz + A), 
R-2 = [V2W1 + xiw 2 + zi(w 3 + X)]w 2 (w 3 - A) + 

+ X 2 ZiW\ + X\Z 2 W 2 + z 1 z 2 (w 3 + A). 

Proof. We need to prove that 

£umuD = €. (3.41) 

It follows from Propositions Q] and [3] that £ C £oU£i. Indeed on £ we have dK = 0, 
dG = ±ab dH . Note also that the system of relations (|3.38p satisfies the conditions of Lemma[5] 
with n — 3, k = 2. Therefore, £ is a smooth two-dimensional manifold with the induced 
Hamiltonian system with one degree of freedom. 

According to Proposition [3] we have D* C £ 2 . Hence, D = D*U£c£. 

Consider the set 91 defined by Eqs. (|3.39|) . First, investigate the cases when these equations 
cannot be solved with respect to yi,y 2 - Suppose that 

UHW2 + Xw 3 = 0. (3.42) 

Then, after several differentiations in virtue of (|3.19p , we come to Eqs. (|3.30p - (|3.33p with 
a = —A 2 . The corresponding points belong to £1. Let 

(W2X1 + Xzi)(wiX2 + Xz 2 ) = 0. (3.43) 

Then the same procedure leads to the system of equations having the only solutions of the form 
(pPg]) . i.e., to the set £. Denote 91* = 9T\(£ U £1). On this set from 1335)) we obtain 

Vi = 7 — : 77 — : — rr[w2(wl + x 1 )(x 2 z 1 w 1 + x 1 z 2 w 2 - 

(W1W2 + Xw 3 )(W2Xi + Xzi)\ 

— X1X2W3 + 2ziZ 2 \) + X2{wiW 3 + Zi)(wiZi — XiW 3 )\— 

-(iii»3 — 2z\Wiw 3 — z 2 )z 2 A 2 ], 
V2 = 7 — : 77 — : — rr[wi(wl + x 2 )(x 2 z 1 w 1 + x 1 z 2 w 2 - 

(W1W2 + AW3)(WiX2 + Xz 2 )X 

-X1X2W3 + 2ziz 2 \) + xi(w 2 w 3 + z 2 )(w2Z2 - a; 2 W3)A- 

— (x2U)l — 2z 2 W2W 3 — z 2 )z\X 2 ]. 

The derivatives of Fi and F2 in virtue of (|3.19p vanish identically after the substitution of the 
expressions Q3.44p . This fact proves that 91* is an invariant set. The Poisson bracket {Fi,F 2 } 
with (|3.44p takes the form 



(3.44) 



{F 1: F 2 } = V2X( Wl w 2 + Aw 3 ) 3/ V(w 2 xi + \zi)(wix 2 + Xz 2 ) C. (3.45) 

Here 

C = i(8s 3 A 2 - r 4 ) J2s 2 - {2h + X 2 )s + p 2 (3.46) 
s 

depends on the energy constant h and the constant s of the partial integral 

s _ xix 2 w 3 - x 2 ziwi - xiz 2 w 2 - Agi£2 ^ 
2A(wiw 2 + Auj 3 ) 

This integral is similar to (|3.37|) and independent of H almost everywhere on 91*. 

The cases when the multipliers (|3.42p or (|3.43|) in (|3.45|) turn to zero are already studied 
above. Zeros of the function (|3.46p have codimension 1. Hence, (I3.45P is non-zero almost 
everywhere on 91*. It follows from Lemma [5] that 91* C £ 2 . Thus, £ U 91 U D C £. To prove 
the equality (|3~4Tj) . we must show that £ C £ U 91 U £>. 
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The points of the set Co described in Proposition Q] satisfy (13.38| . According to Proposi- 
tion [5] the set £1 can be represented as £n U £12, where £n consists of the trajectories (|2.14p 
and £12 is defined by the system (|3"7M)) - (j3T2"?|) . On the trajectories (|2?T3|) we have (EHgj) . It 
is easily checked that the points given by Eqs. (|3.24[) — (|3.26|) under the condition (|3.27[) satisfy 
both systems ([335)) and (j3~¥)]) . Therefore, £ U <£ u C £ and £ 12 C 9t n D. 

Consider now the set £2- To investigate the dependence of G, H, K, introduce the function 
with Lagrange's multipliers. It follows from Propositions Q] and O that on £ 2 the differentials 
dK and dH are linearly independent. Then the multiplier at the function G is always non- 
zero and can be chosen equal to any non-zero constant. It is convenient to take the function 
2G + SK + (T — p 2 )H, where S and T are Lagrange's undefined multipliers. The condition 

2dG + SdK + {T~p 2 )dH = (3.48) 

is preserved by the phase flow. Applying the Lie derivative we obtain 

SdK + fdH = 0. 

Since rankjdG, dK, dH} — 2, this linear combination of the differentials is proportional to the 
left part of (|3.48|) . It means that, at the points of £2, 

5 = 0, r = o. 

Thus, the functions S and T are the partial integrals on the submanifold £2. According to 
Lemma H] rewrite Eq. (|3.48p as the system 



£2(2/2 + 2S)wi + 2S(w\w 2 + Xw 3 )w 2 + 

+ (T - z x z 2 - 2SX 2 )w 2 + x 2 z x w 3 + (y 2 + 2S)z 2 X = 0, 
xi{yi + 2S)w 2 + 2S{w\w 2 + Xw 3 )w x + 

+ (T - z x z 2 - 25A 2 )wi + x x z 2 w 3 + (y x + 25)ziA = 0, 



(3.49) 



(T - x x x 2 )w 3 + x 2 z 1 wi + x 1 z 2 w 2 + (2Sw x w 2 + ziz 2 )X = 0, 

Tz\ + x x z 2 w\ + {(xix 2 — 2z 1 z 2 )wi + (y x z x + x x z 2 )X + x x y x w 2 ]w 3 — 
— (y x z x + xiz 2 )wiw 2 + xi(yi + 2S)w 2 X- 

-[x 2 zi + {y 2 + 2S)z 2 ]w\ + [(y 2 + 2S)y x - 2z 1 z 2 ]w 1 X + y x z x X 2 ~ 
-[{y 2 + 2S)x x + z 2 ]z 2 =0, 

Tz 2 + x 2 z x w\ + [{xix 2 - 2z x z 2 )w 2 + (y 2 z 2 + x 2 zi)X + x 2 y 2 w x ]w 3 - 

-(U2Z 2 + x 2 z x )w x w 2 + x 2 (y 2 + 2S)w x X— (3.50) 

~[x x z 2 + (y x + 2S)z x ]w 2 + 

+[(y x + 2S)y 2 - 2ziz 2 ]w 2 X + y 2 z 2 X 2 - 

-[(y x + 2S)x 2 + z 2 ]z x =0, 

(T - xix 2 )(yi - 2/2) + 2(y 2 + S)x 2 w\ - 2{y x + S)x x wj + 
+2(x 2 ziwi - xiz 2 w 2 )w 3 + x 2 z\ - x\z\+ 
+2(y 2 z 2 w x - yiZiw 2 )X = 0. 

It follows from Proposition [3] that this system is valid at the points of the set D*. To find 
all other cases suppose 

UiU 2 ? (3.51) 

and express y x , y 2 from (|3.36p : 

2/1 = ~ t ^i 2U 2 ~ [wiw 2 (w 3 - A) + W 2 Z 1 - W X Z 2 \U X - 

2w x w 2 {w 3 — X) 

- 2[wiz 2 (w 3 - A) 2 + {x 2 w\ + z ± z 2 + 2Xwiz 2 )(w 3 - X) + 
+ x 2 z x w x + x x z 2 w 2 + 2Xz x z 2 }}, 



2/2 = T rr{ 2 ^2 + [w\w 2 {w 3 - A) + w x z 2 - w 2 z 1 }U 1 - 

2w x w 2 (w 3 — X) 

- 2[w 2 zi(w 3 - A) 2 + (x x w 2 + z\z 2 + 2Xw 2 z\){w 3 - A) + 
+ x 2 z 1 w 1 + x\z 2 w 2 + 2Aziz 2 ]}. 



(3.52) 
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The determinant of the system p. 491) with respect to T, 25 is equal to 

A = X±W 2 — X 2 w\ — {z 2 W\ — ZiW 2 )X. 

If we suppose that A = on some time interval, then the sequence of the derivatives of this 
identity in virtue of (|3 . 1 9f> leads to (|3.38p . i.e., to the points of JjUCi. Consider then 



A ^ (3.53) 



and find from Eqs. (|3.49|) 



5 = — [x 2 y2w\ - x x yiwl + (x 2 z 1 wi - x 1 z 2 w 2 )w 3 + 

+ (y 2 z 2 W! - y 1 z 1 w 2 )X], (3.54) 

T=^[A 1 B 1 -A 2 B 2 ]. 



Here 



Ai = {x 1 w 2 + Xzi)y 1 + (xiw 3 - ziwi)z 2 , 
B\ = {u>2 + x 2 )wi + Xw 2 (w 3 — A) + Xz 2 , 
A 2 = (x 2 wi + \z 2 )y 2 + {x 2 w 3 - z 2 w 2 )zi, 
B 2 = [w\ + xi)w 2 + Awi(w 3 — A) + Azi. 

Substitute (|3 . 52[1 and (|3.54jl into (|3.50p to obtain the system of four equations of the type 
Ei = [i = 1, 4), where E^ = ai 2 U\ + anL \ +<iioU 2 with some polynomials fly. For the proof 
of the theorem, there is no need to use all equations of this system. It is enough to consider, 
for example, the zero points of the resultant of two simplest functions E\,E± with respect to 
U 2 . We obtain Awiw^iRA = 0, where 

R = \wiw 2 (xiw 2 + \zi)(x 2 wi + \z 2 )Ui - {wiw 2 [x 2 ziwi + XlZ 2 W 2 - 

- xix 2 {w 3 - A) + 2Aziz 2 ] + (ziz 2 w 3 + x 2 z\W\ + xiz 2 w 2 )\ 2 + ziz 2 \ 3 } A. 

Due to (|3.51[) . U\ 7^ 0. At the points of C\£ the product w\w 2 is not identically zero. The 
set £\(£UD) is, obviously, preserved by the phase flow. Therefore, (|3.53[) and (|3.19|) imply 
R = 0, R' = 0. These equations are linear in yi,y 2 and yield the expressions (|3.44p satisfying 
Eqs. (13391) . Thus, £\(£ U D) C STl. □ 

Remark 4. The system (f3"H))) follows from Eqs. (I33B1) and ([335)) . In particular, D* = D\£. 
Then from Proposition [3] and Lemma [5] we obtain that this manifold lies completely in £ 2 . 

Remark 5. The integral 5 given by J3.54[) . when restricted to the set 0^ has the form (|3.47[) . 
Indeed, it is enough to substitute (|3.44|) into the first formula (I3.54|) to obtain (|3.47[) . On the set 
D the same function 5 in the substitution of (|3.35[) and (|3.52p takes the form (|3.37[) . Therefore, 
the use of the same notation in (|3.37[) and (|3.47[) is correct. The expressions for T can also 
be simplified. On the set 01 we have T = 2 A 2 5, i.e., this function does not give rise to a new 
partial integral independent of 5. On the contrary, at the points of the set D we have 

T = x\x 2 + z\z 2 — 2w\w 2 S. 

In the case A = the same expression with the corresponding function 5 is the partial integral 
independent of 5. The equations of the integral manifold defined by the pair 5, T lead to the 
separation of variables on O [M] . 

4 The bifurcation diagram 

The Lax representation for the considered problem found in [2] can be written in the form 

L' = LM - ML, (4.55) 
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where 



L = 



1 \ 


■I 2 


— 2llU 


z 2 




w 3 


n 


il'2 














X 


X 




2 




~2 


Xl 


-2A 


Zl 


2wi 







w 3 





X 

-2wi 


Z-2 


X 

-2w 3 


-£-4* 


, M = 


Wi 


2 







X 


X 




2 




T 




2w 2 


h 4x 


2w 3 







w 2 


— x 


X 




X 








2 







Wl 

' 2 
x 

W3 

' 2 



Here x stands for the spectral parameter, the derivative in (|4. 55[) is calculated in virtue of the 
system p,19j) . The equation for the eigenvalues n of the matrix L defines the algebraic curve 
associated with this representation [8]. Let s = 2x 2 and let h, k,g be the arbitrary constants 
of the integrals (|3.22|) . The equation of the algebraic curve takes the form 



2rP_ _ {2h + A 2) + 2s] + 4[ !_ + ? (4 _ 2p 2 h _ p 2 x 2 ) + 

s s s 

+4{k + 2\ 2 h) -8A 2 s] = 0. 



(4.56) 



It is natural to suppose that the bifurcation diagram of the momentum map (|3.16p is 
included in the set of values (g,k,h) such that the curve (|4.56p either have singular points or 
is reducible, i.e., the left part of Eq. (|4.56|) splits into the product of some rational non-trivial 
expressions. In this way we can guess the result of the following statement. Nevertheless, 
to obtain the complete proof of it, we must fulfil the calculations on the above found critical 
manifolds. 

Theorem 3. The bifurcation diagram of the momentum map G x K x H is included in the 
union of the following (intersecting) subsets ofM. 3 (g,k,h): 
1 ) the pair of straight lines 



+ ■ 



2) the surface 



3) the surface 



= {a + bf, 

X 2 

= -ab(h- —); 



AX 2 s-2X 2 h- 



r 

4^ 



= (a-b) 2 , 
X 2 

= ab{h- —); 



r 

4? 



s G 



(4.57) 



(4.58) 



k = '3s 2 -A(h- 



2 
A 2 



(h 



4s 



S G 



4s 2 
^\{0}. 



(4.59) 



Proof. Let ( G £. Substitution of the values z\ — Zi = into (|3.20[) and (|3.2ip yields x±x 2 — 
(a ± b) 2 , yiy 2 = (a T b) 2 . Then from (|3.22p . (|3.38p we obtain the equations defining the lines 

Let C G Take the constant of the partial integral (|3.47p for the parameter s in (|4.58p . 

substitute the expressions (|3.22p for the corresponding constants, and fulfil the change (|3.44p . 
Then both Eqs. (|4"35)) become the identities. Therefore, J(9t\£) C Ti. 

The inclusion J(D\£) C T2 is proved in a similar way. We take the constant of the 
partial integral (|3.37p for the parameter s in (|4.59p and fulfil the substitution (|3 . 52[) with 
U± = U 2 = 0. □ 
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Remark 6. Note that the shift of the energy level h = h — A 2 /2 makes the equations of the 
lines r± and the surface T2 independent of A. Thereby obtained equations are identical with 
the corresponding equations of the case A = [H]. The surface Ti is obtained as a perturbation 
(with respect to A) of two tangent to each other sheets of the bifurcation diagram of the case 
A = 0, i.e., the plane k = and the slanted parabolic cylinder (p 2 h — 2g) 2 — r A k = 0. Thus, 
it is easy to view the evolution of the Appelrot classes [50] of the S. Kowalevski case in the 
process of two-way generalizations — adding the second force field and, afterwards, the non-zero 
gyrostatic momentum. 

The equations given in Theorem [3] are in the following sense convenient. Let us fix the 
energy constant h. Then we obtain the parametric equations of a one-dimensional set in the 
plane (g, k) (with the finite number of singular points). This set is the bifurcation diagram E^ of 
the restriction of the pair of integrals G, K onto the iso-energetic surface {H = h} C P 6 , which 
is always compact. In particular, all diagrams E^ lie in the restricted area of the (<?, A:)-plane and 
are easily drawn numerically. The analytical investigation of the types of the diagrams E/j with 
respect to the essential parameters (b/a, X/y/a, h/a) is a necessary but technically complicated 
problem. Nevertheless, it must be solvable. Indeed, the set of double points and cusps of 
the curves IY2 in the (g, fc)-plane is easily defined and investigated analytically. Moreover, the 
values of the first integrals on the periodic motions (13.24j) - ()3.27p define the points of transversal 
intersections Ti n IV This fact, at least, guarantees that the numerical algorithm can be built 
for effective calculation of knots of one-dimensional cell complex E/j for any h. In turn, it 
should be possible to find all cases of bifurcations of the set of these knots with respect to the 
parameters defining the above set of periodic motions. 
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